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Abstract 

A generalized filter construction is used to build an example of a non- 
MRA normalized tight frame wavelet for dilation by 2 in L 2 (R). This 
example has the same multiplicity function as the Journe wavelet, yet has 
a C°° Fourier transform and can be made to be C r for any fixed postive 
integer r. 
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1 Introduction 

A central problem in the history of wavelet theory has been the search for smooth 
well-localized functions whose dilates and translates can be used to obtain all 

*The first two named authors were supported by a US-NSF Focused Research Group (FRG) 
grant. 
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members of a specified function space. We will use here the simplest definition 
of an orthonormal wavelet, by taking our context to be dilation by 2 in L 2 (M.): 

Definition 1. tp € L 2 (M.) is an orthonormal wavelet «/{V'j,fe = V2 J ^(2^ x — 
k)}j,kez form an orthonormal basis for L 2 (R). 

Many of the techniques and results used in this paper apply in much wider 
contexts, to dilation by an expansive matrix in L 2 (R n ), and to allow for multiple 
wavelets. However, for the purposes of creating the example that is the main 
thrust of this paper, this more restrictive context will suffice. 

The earliest example of an orthonormal wavelet was given by A. Haar in 
1910 PH]- Here the wavelet is defined in terms of characteristic functions of 
bounded sets, 

and thus is well-localized but not at all smooth. Another famous classical ex- 
ample, the Shannon wavelet, has Fourier transform of the wavelet equal to a 
characteristic function of a bounded set, 

^=X[- 1) -i)u[i,i)> 

and thus has the opposite problem of being smoooth but not well-localized. 
In 1988, Ingrid Daubechies ^7] gave a construction technique for orthonormal 
wavelets for dilation by 2 in L 2 (R) that are both smooth and well- localized. She 
described how to build, for any integer r > , a wavelet that is C r and has 
compact support. It can be shown that no wavelet can be both C°° and have 
compact support (see e.g. [201 ) - 

Daubechies' construction technique used filters in a multiresolution structure 
called an MRA. Most classical wavelets are associated with such a structure. 
J.L. Journe gave the first known example of a wavelet that can not have an 
associated MRA (see [22]). Journe's wavelet, like the Shannon wavelet, is a 
wavelet set wavelet. That is, its Fourier transform is a characteristic function. 
Since the announcement of the Journe wavelet, many other non-MRA wavelets 
have appeared in the literature (see e.g. ^5]> HI most of them also wavelet 
set wavelets. These wavelets are smooth but not well-localized. In fact, it can 
be shown that every orthonormal wavelet ip such that \ip\ is continuous and 
\i>((i)\ = 0(|Cr^" e ) at oo for some e > must be an MRA wavelet (see [20]). 
Thus, in the quest to find wavelets that are both well-localized and smooth, 
non-MRA wavelets would seem to be not very useful. However, if we expand 
our definition of wavelet slightly, the situation changes. 

Definition 2. The set of functions {r/j} is a frame for L 2 (R) if th ere exist con- 
stants A and B such that for each f £ L 2 (R) we have A||/|| 2 < \{f\rij)\ 2 < 
i?||/|| 2 . The set {r]j} is a normalized tight frame or Parseval frame if A = B = 
1. 

Definition 3. ip S L 2 (R) is a (normalized tight) frame wavelet or a Parseval 
frame wavelet for dilation by 2 if {4>j,i = \[ 7 H ^(f2? x — I)} form a normalized 
tight frame for L 2 (R). 
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Note that a normalized tight frame can exhibit redundancy, and therefore 
need not be a basis. Indeed, it can include as one of its elements. However, 
a normalized tight frame {r)j} does have the property that every / £ L can 
be recaptured from its coefficients, / = Yl{f>Vj}Vj- (See, e.g. |2J]|). Thus 
replacing orthonormal wavelets by frame wavelets does not give up the essential 
property of determining all functions in L 2 in terms of translates and dilates of 
the wavelet. 

Our main goal in this paper is to produce an example of a non-MRA normal- 
ized tight frame wavelet for dilation by 2 in L 2 (R) which is both well-localized 
and smooth, in the sense that i/j £ C r (for any fixed positive integer r) and 
ip £ C°° . The existence of such a wavelet is somewhat surprising in view of 
the result mentioned above for orthonormal wavelets. While several researchers 
have produced examples of frame wavelets that exhibit smoothness and good 
localization properties (See e.g. |T]), the example in this paper stands out in 
its relationship to a non-MRA multiresolution structure with a known multi- 
plicity function. For such a non-MRA wavelet, the multiresolution structure 
will be replaced by a generalized multi-resolution structure called a GMRA, 
defined in which we review in Section 2 of this paper. We note that similar 
generalizations of multi-resolution structures have also been defined in [H] and 

EH- 

The Fourier transform is an essential tool in analyzing wavelets; we will use 
it and a generalization in the GMRA context to define filters used in wavelet 
construction. We take as the definition of the Fourier transform, 

f(x) = [ f{t)e- 2 ™ l dt. 

JR 

Note that on L 2 (R), translation by n becomes multiplication by e 2mnx , and 
dilation becomes Sf(x) — ^j/(§ )■ We review the filter construction process 
and its generalization in Section 3 of this paper. In Section 4 we build our main 
example. 

2 Generalized multi-resolution analyses 

As mentioned in the introduction, most classical wavelets have an associated 
multi-resolution structure of the following type, first defined by S. Mallat |22j : 

Definition 4. A Multiresolution Analysis (MRA) in L 2 (R) is a collection of 
closed subspaces Vj that have the following properties: 

1. V, C V j+1 

2. V j+1 = {6(f) = V2f(2x)} feV] 

3. L)Vj is dense in L 2 (R) and C\Vj = {0} 

4- Vo has a scaling function (j) whose translates form an orthonormal basis 
for Vo 
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However, the Journe wavelet (see \'22\). 



1> = X[-^,-2)U[-i -| 



f)U[|,i]U[2 M). 



has no associated MRA. This can be shown using the dimension function, 



first introduced by P. Auscher 0. Because an MRA requires the existence of a 
scaling function, it turns out that a wavelet has an associated MRA if and only if 
the dimension function is identically equal to 1 (see e.g. PHI)- By removing the 
scaling function condition, we are able to associate a multiresolution structure 
with the Journe wavelet. More generally, given any orthonormal wavelet ip € 
L 2 (R), the subspaces {Vj} defined by Vj = the closed linear span of {tpi,k}l<j, 
do determine a generalized multiresolution structure, according to the definition 
below 0: 

Definition 5. A Generalized Multiresolution Analysis (GMRA) is a collection 
of closed subspaces {Vj}j^z of L 2 (R) such that: 



4- Vo is invariant under translation. 

The definitions of MRA and GMRA differ only in condition (4): An MRA 
requires that Vo has a scaling function <j> such that translates of <p form an 
orthonormal basis for Vo, while a GMRA requires only that Vb be invariant 
under translation by the integers. In spite of this difference, it is shown in 
jS] that a GMRA has almost as much structure as an MRA. Translation is 
a unitary representation of Z on Vq, and thus is completely determined by a 
multiplicity function to : [-|,|) i~> {0, 1, 2, • • • , oo} describing how many times 
each character occurs as a subrepresentation. The multiplicity function has 
been shown |25j to be equal to the dimension function defined in Eq. l|T]l. By 
writing V\ = Vo ® Wo, representation theory can be used (see 0) to show that 
the GMRA has an associated orthonormal wavelet if and only if the multiplicity 
function satisfies a consistency equation: 



A GMRA is an MRA iff to = 1, and Journe's famous non-MRA wavelet example 



oo 




(1) 



fcez i=i 



1. V 3 c V j+1 

2. V J+1 = {6(f) ee V2f(2x)} feVj 

3. UVj dense in L 2 (R) and nVj = {0} 




has 
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In every GMRA, there is a unitary equivalence between translation on Vo and 
multiplication by exponentials on ®L 2 (Sj), where Sj — {x : m(x) > j}. This 
unitary equivalence plays a role here similar to that of the Fourier transform in 
the classical MRA case. It ensures that a GMRA whose multiplicity function 
is finite a.e. has generalized scaling functions, (pi, (p2, • • • such that {(f>i(x — I)} 
form a normalized tight frame for Vo- If m is bounded, with maximum value 
c, then c generalized scaling functions are required. If the GMRA is associated 
with a wavelet set wavelet, these functions can be chosen of the form <pi = XEi, 
where XE f is the characteristic function of the set E it chosen to be congruent 
mod 1 to Si, and so that LiEi C 2 (UEi). (For more details see Theorem 3.5 
of |7].) For example, the Journe wavelet has generalized scaling functions (pi = 

X[^,l)u ± [i,i) and <fo = x ±[ i,f) 



3 Filter constructions 

We begin by recalling the classical filter construction developed by Mallat [22 
and Meyer j2H], and used by Daubcchies ^B] to build C r wavelets with compact 
support. Thus, suppose we have a single wavelet ip for dilation by 2 in L 2 (M), 
with an associated MRA, and so a scaling function <f> whose translates form an 
orthonormal basis for Vo ■ 

Because <f> and tb are in Vi, we can write <f> and ip in terms of exponentials 
times the dilate of <f>. That is, there must exist periodic functions (with period 
1) h and g such that 

fa) - ^(f)£(f ) (2) 

and 

to = -^= 5 (| m x - ). (3) 

The functions h and g are called low and high pass filters. The roots of 
this name can be seen by noting that for the Shannon wavelet, h — \/2X[_ i i) 

and g = \/2X[-i,-i)u[i,i) d° indeed act by filtering out all but low (for h) 
or high (for g) frequencies. Because of the orthonormality conditions satisfied 
by translates of <j> an d ^ a U filters defined by Eqs.J2J) and (J3J) must satisfy 
orthonormality-like conditions: 

\h{x)\ 2 + \h{x + \)\ 2 = 2 (4) 



I^)| 2 + l5(^+2)| 2 = 2 



and 



h(x)g(x) + h(x ■ 



1 ^ < 1 t 



(5) 
(6) 
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The classical filter techniques for building wavelets reverse this process of 
finding filters from wavelets. First note that if h is any periodic function that 
satisfies Eq.(QJ, the function 



g(x)=e 2 ™h(x + -), (7) 

and h together satisfy all three orthonormality conditions an d ©■ (Other 

choices for g are possible.) The following classical theorem gives conditions 
under which the Fourier transform of a scaling function (p can then be built by 
iterating Eq.0, so that the Fourier transform of a wavelet, ip, can be defined 
using Eq.(J3J). 

Theorem 1. Let h and g be C 1 functions that satisfy the orthonormality condi- 
tions 0), and 0). Suppose, in addition, that h is nonvanishing on [—\, \), 
and |/i(0)| = y/2. Then: 

fa) = n -^v-**) 

3 = 1 VZ 

is a scaling function for an MRA, and 

fa) = ±4 )?(f). 

is an orthonormal wavelet. 

A good introductory description of these constructions appears in [21] • 
The theorem's requirements that h satisfy \h(0)\ = y2 and that h be in 
C 1 are natural restrictions in order to make the infinite product converge. The 
condition that h is nonvanishing on [— \, \] appears less natural; it is used 
in the proof to ensure L 2 convergence of the infinite product and thus the 
orthonormality of the translates of (p. A famous example due to A. Cohen 
12 showed that this nonvanishing condition cannot be entirely removed, as its 
removal can lead to functions <fi and ip whose translates are not orthonormal. 
Cohen took h = 1+6 ^ — , which resulted in a stretched out version of the Haar 

v2 

scaling function and wavelet, (f> = gX[o,3) an d ip — \{X\o,&) ~ X[8,3))- Cohen's 
tp is not a wavelet, since its translates are not orthonormal. However, the 
classical Theorem^can be extended to accommodate this and similar examples 
if we allow (normalized tight) frame wavelets. The following generalization of 
Theorem^ was first proven by Lawton in |21|. with a generalized form appearing 
later in Bratteli and Jorgensen's book 

Theorem 2. Suppose h, g are periodic Lipschitz continuous function in L 2 (R), 
which satisfy \h(0)\ — \2 and the filter equations 0) and |^). Then the 
construction (p(x) = Y]!jLi ^^( 2_ '' ' x ) P r °duces an L 2 function cp (whose trans- 
lates are not necessarily orthogonal), and the function ip(x) = ^fl'(§) < ^(§) * s a 
(normalized tight) frame wavelet for dilation by 2 in L 2 (M). 
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The classical filter construction that was used in both the original theorem 
and this extension was generalized to GMRA's in 0|. Just as in the classical 
case, it is natural to begin by building filters from wavelets, and then determine 
conditions under which the process can be reversed. Accordingly, suppose first 
that ip is a non-MRA orthonormal wavelet for dilation by 2 in L 2 (R) whose 
multiplicity function is bounded by c. We have seen in Section 2 that there 
then exist generalized scaling functions <pi, ■ ■ ■ , </> c . Since <fi and ip are in V\, we 
can show (see U) there exist periodic functions hij and gj, supported on the 
periodization of Sj , such that 

&(z) = ^£Mf)^(f) (8) 

3=1 

and 

^) = 4l>(f)^(f)- (9) 

3=1 

These generalized filters gj and ftjj satisfy orthonormality-like conditions that 
are generalizations of the classical conditions , © and JJjJ : 



3 = 1 

EMf)| 2 + b(^)| 2 = 2, (11) 

3=1 

and 

EMfta4) + ftij(^)fc(4^) = °> (12) 

3=1 

For example, for the Journe wavelet these filters are defined (see by : 
hi,i = V2x[_|,_j) U (-^^)u[|,f ) 
Zn,a = 0, 

^2,2 = 

ai = 44,-i)u[H) 

.92 = V2X[-i,i)- 
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To use generalized niters to build new wavelets, we reverse this procedure, 
just as in the classical case. In order to first build filters, we use functions on 
the disjoint union of the Sj's whose values are v2 times unitary matrices, with 
different dimensions for different values of x. We need the values of the filters to 
be \[2 times unitary matrices in order to satisfy the generalized orthonormality 
conditions (|l(Jfl . Q11JI- and i|12|l . The matrices of filter values have different 
dimensions depending on how many of the sets Sj the point x and its preimages 
are in. Once we have the filters, we build the generalized scaling function using 
an infinite product of matrices that comes from the iteration of Eq.JHJ). The 
wavelet is then produced by Eq.Q. Conditions that make this possible are 
described in the following generalization (see 0) of the Lawton theorem : 

Theorem 3. Given a multiplicity function m for a GMRA, suppose {hij} 
and {gj} are periodic functions that are supported on the periodization of Sj = 
{x : m(x) > j}, Lipschitz continuous in a neighborhood of the origin, and 
that satisfy the three generalized orthonormality conditions ll(J\) . and \l 6 3jt . 
Suppose in addition the hi_j satisfy the generalized lowpass condition \hi j(0)\ = 
v / 25(i ! i)(5(j.i) • Write H for the matrix (hi,j). Then the components of 
YYf°—i ~^/2-H(2~ k x) converge pointwise to 1? functions. If we let {4>i\1 = i be the 
first column of this product, then the translates of {4>i} determine the core sub- 
space Vo of a GMRA, and 

^)^-L ffl( |)^(|) + 92 (|)^(|) 

is the Fourier transform of a normalized tight frame wavelet on L 2 (M) . 

The proof, like that of proceeds by using matrices of values of the 
niters to define partial isometries that satisfy relations similar to those defining 
a Cuntz algebra as in |14|. We will use this theorem in the next section to build 
our main example. 



4 The Example 

For any fixed positive integer r, we will now use Theorem to build an example 
of a C r frame wavelet with C°° Fourier transform, whose multiplicity function 
is that of Journe: 

r 2 I6 [-U) 

m(x) = i 1 xe±[±,f)U±[f,±) 
[ otherwise 

To do this, we will first choose low-pass filters that force the functions 4>\ and </>2 
in the infinite matrix product to both be smooth and have rapid decay. We will 
then define the high-pass filters, and show that the resulting wavelet inherits 
the properties we seek. 



8 



From the equation for to, we see that the sets Sj — {x : m(x) > j} used in 
the filter construction are here given by S\ = [— |, |)U±[|, i) and S2 = [— j, j). 
We will define the filter functions hi_j, 1 < i,j < 2 in terms of a classical filter 
p. Accordingly, with r fixed as above, let p be a periodic (with period 1), 
real-valued function that satisfies: 

1. peC°°(R) 

2. p(x) = p(-x) 

3. \p(x)\ 2 + \p(x+±)\ 2 = 2 



4. pW(x) = Oforxe [i,^]U[f,i], fc>0. 

5. |p( r+2 )(a;)| < 1 when \x - 



j|| < or |x- 



±1 < — 

7 I v 112 



To see that such a function exists, note that it suffices to define p e C°° ( , j] U 
tM' iD' satisfying conditions (4) and (5), and with p(^g) — p{\) — p{%) = 1 
and p {k) {-^) = = P^H^) = for k > 1. We can then define p on 

[0, jg] U [\, Ig] by p(ir) = ^2 - - .t) 2 . We then extend p to [-§,0] using 
the requirement that p be symmetric, and finally extend p to the whole line 
using periodicity. Figure 1 below shows a C 1 approximation to p. 

V2h~ 



1 

V2 



Figure 1: AC 1 approximation to p 
We now use this function p to build the low-pass filters {hij}: 

_j P (x) ze[-f,f) 



/ii,i(a;) 







otherwise 



(13) 
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hi. 



l 2,l 



'p{x+\) x£[-\,\) 
otherwise 

_ fy/2 a;€±[?,£) 







7' 2' 

otherwise 



(14) 
(15) 



& 2 ,2(a0 = (16) 
Figures 2 through 4 show C 1 approximations to the three nonzero low-pass 



filters. 
V2h 



l 



Figure 2: An approximation to h\,\ 

It is easy to check that these functions satisfy the filter equation (|10|) . For 
example, to check that 

1, 



iKibr + i^c^r + |/n,i(-5-)r + i^ 2 (^)i 2 = 2 X[ -, f] 



'2" 1 i "^2" 1 ' v 2 
we note that the left hand side gives 

|2^|p(x±l)|2 + 

+ 6 + 

o + b(^)| 2 



u±[ 



b(f)l 2 
b(f)l 2 



2forze 







for x e ±[f , 7 
2 for x e ±[f,i 



2 ,i 2 

The other cases of (llOf) follow similarly from the definitions of the ftjj and the 
properties of p. It is also easy to check that the {hij} have proper support and 
satisfy the low-pass condition, hi,i(0) = v2 and /ii,j(0) = for (i, j) ^ (0,0). 

As in Theorem^ we will define the functions <j>i and <j>2 (and thus the GMRA 
that they determine) in terms of the low-pass filter {hij}. First we will need 
the following technical lemma: 
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V2 



1 

V2 



Figure 3: An approximation to 



Lemma 1. For every integer n > 1, the (3n) partial product of the filter matri 
whose nonzero components have the following support: 



support(a n ) = 



(s»Z + (hf , f ] U ±[|, |] (J ^ 2 A k (J UZ^B* 
(J(4-8"- 1 Z + yl„_ 1 ) 



support(c n ) = 8™Z + U^Cfc 



where 



Bh 



support{d n ) = 8™Z + C n -i 



{ a ,}e{i,-ip+i j=0 

U 

{a J }e{l,-l}*+ 1 J=0 



and C fe = |J ^ aj8 J ' + [- 
{ 0j }e{i,-i}*+i j=o 



1 ii 
7'7 J 
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V2\ 



1 

V2 ' 



Figure 4: An approximation to /i2,i 



Proof. The proof is by induction on n. Both the n = 1 case and the induction 
step follow directly from the observation that 

AAi,i(#) h h2 (§)\ _ (h^x) \ 
/J V>2,i(#) J \t^(x) t 2 , 2 (x)J 

where 

= ^i,i(|)^i,i(|)^i,i(|)+/ii, 2 (|)^2,i(|)^i,i(|)+^i,i(f )^i,2(|)^2,i(|) 
has support 

8Z + {[-I, 2 -] U ±[^, i]) |J 4Z + 2 ± [1 H] |J 8Z + 4 ± [| 



and both 
and 

have support 



t 2 ,i(x) = h 2A (-)h 1A (-)h 1A (-) 



XXX 

t 2 , 2 {x) = h 2A (-)h 1A {-)h h2 (-) 



8Z± [-,-]. 
L 7' T 



□ 
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Remark 1. Since J3y shows that the infinite product of filter matrices, 




/2 ) 



) 



converges pointwise, Lemma^\ gives the support of the elements of this matrix. 
In particular, the support of the upper left element (which will be used to define 
fa) is [-f , f] U ±[f , f] IJUfco^fc U U fco B fc' a7ld </le support of the lower left 
element (which will be used to define fa) is Uj? =0 Ck- Note that these supports 
are disjoint. 

We will use this result to characterize the GMRA for the wavelet we will 
build. 

Proposition 1. The functions fa and fa given by the upper left and lower 



termine a GMRA with multiplicity function the same as the Journe wavelet. 

Proof. Define Vo to be the closed linear span of exponentials times fa and fa , or 
in other words, define Vo to be the closed linear span of translates of fa and fa. 
A GMRA can be built from this core subspace Vo in a natural way, by letting 
Vj = = f(Vx)} fev - By 5 , the resulting structure does satisfy 

all the components of Definition [5] To find the multiplicity function associated 
with this GMRA, we note first that LemmaHshows that the supports of fa and 
fa are disjoint, so that the multiplicity function is the sum of the multiplicity 
functions determined by the two functions separately. By y], the multiplicity 
function associated with fa is the characteristic function of the support of the 
periodic function Per fa {x) — J2i<ez \4>j{ x + 0| 2 - From Lemma Q we see that 
these sets are indeed the Journe sets Si for fa and 5*2 for fa. □ 

Remark 2. Neither Per fa nor Per fa is bounded away from on its support, 
so that by ]9], the translates of fa and fa do not even form a frame for their 
closed linear span Vq. (This is true in the Cohen example as well.) 

Figures [3 and below show C 1 approximations to the graphs of fa and fa 
near the origin. As shown in Lemma^ both functions have unbounded support, 
but the portions of these functions that do not appear in Figure [5] and Figure [6] 
have values less than .01 and .002 respectively. We will use fa and fa to define 
our C r frame wavelet with C°° Fourier transform. To do this we will need the 
following properties of fa and fa: 

Proposition 2. The functions fa and fa defined in Proposition^ are C°° . 

Proof. We fix an a and show that fa and fa are C°° in a neighborhood of a. 
Because /i.2,1 is on (— all but a finite number of the lower triangular 
products of three successive factors described in Lemma ^ are in fact diagonal 



left elements of the infinite product 
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with single terms on the diagonal. Thus, each of the entries in the first column 
of the infinite product matrix is the sum of only a finite number of terms. Since 
= \pl in a neighborhood of 0, each of the terms has only a finite number of 
factors not equal to 1. It will suffice to show that each of these finite products 
is C°°. 

By construction, we have that the hij are all C°° everywhere except for 
hi,i at n ± |, h% t 2 at n ± h, and h 2 ,i at n ± | for n £ Z. We will show that 
whenever one of these discontinuities occurs as a factor in the infinite product, 
it is cancelled by a following factor that is at the point of discontinuity. From 
the formula for the tij given in Lemma ^ we see that any term in the infinite 
product that contains a factor of h 2 .\(n ± j) must also contain a factor of one 
of the forms hi t i(n ± |), h 1 ,i(n ± ^), hi^(n ± ^), or h\, 2 {n ± |). The last 
three possibilities are in a neighborhood of the point in question, so if we have 
a discontinuous factor of it is either cancelled out by a factor, or we also 
have a factor of h\^{n ± |) with a smaller n. Similarly, any term in the infinite 
product that contains a factor of hx,i(n ± |) must also contain a factor of one 
of the forms h\^{n ± y|), hi^(n ± j), hi^in ± 75), or hx,i{n ± ^). The first 
three of these possibilities are at the points in question, so any discontinuous 
factor of hi t i is either cancelled out by a factor, or is followed by a factor of 
hi^{n± j) with an equal or smaller n. Finally, any term in the infinite product 
that contains a factor of hi^(n ± j) must also contain a factor of either the 
form h 2: i(n ± or the form h 2: i(n, ± |). The first of these possibilities is at 
the point in question; the second possibility throws us back into the first type 
of discontinuity we considered above, but with a smaller n. We can repeat the 
above sequence of arguments if necessary, noting that each succeeding factor is 
evaluated at a point half the distance from the origin as its predecessor, so that 
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.06 



.03 



Figure 6: AC 1 approximation to <f>2 



the chain above must eventually end with a factor of hx,i and thus with 0. 

Thus each of the finite products that occur as terms in <pi and <p2 has the 
form of product f(x) — /i(x)/2(x) in a neighborhood of o, where /i(x) is C°° 
on (a — e, a + e), with its value and all its derivatives, f[ n \a) = 0, n > 0, and 
/2(x) is C°° on (a — e, a) U (a, a + e), with a jump discontinuity at a, but with 

0, n > 1. Thus, 
□ 



all derivatives satisfying lim :E ^ a - /| (x) — . , 
/ is C°° on (a - e, a + e) with / (n) (a) = 0, n > 



lim,^„+ /a (x) 



The next result shows that the <^>j vanish rapidly at oo, to an extent deter- 
mined by the integer r chosen at the beginning of the example. 



Proposition 3. The functions <p\ and <f>2 defined in Proposition^satisfy 
as x — ► oo. 



|r+l 



4>i(x) 



Proof. Using Lemma ^ we see that it will suffice to show that given an e > 
0, there exists N such that |x| r+1 |0i(a;)| < e for x G A n U B n U C n with 
n > N. Recall from that lemma that x € A„ U £?„ U C„ has the form x — 
d (^2^ =0 aj8 J + ±ro^ , where de 1,2, 4, ejj = ±1, and where < ro < j. For any 

such x, 4>i(x) is the sum of less than 3™ +1 terms of the form YiuLi ^^ndiiw)' 
since once x < j, we have ft»2,i(x) = 0, so that the triangular matrix de- 
scribed in Lemma ^ is actually diagonal with single term elements. To see 
how the argument of the j t changes as we move from 
that if y of the form y = J2j=o a k+j& + a k-i r k with j — 5 < r < 
I = 4 a k+j+i& + flfeT-fe+i) , | = 2 1 a k+3+i8 3 + a k r k+1 j , 

| = Yll=o 1 a fe+i+i8"' + flfe^fe+i. The remainder terms, g? ■ rje+i, take one of 



% to we note 

i, then 

and 
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two forms, depending on the value of ak-i- If cik-i = +1, then Ar^+i is within 
| of |, so that | — | < rk+i < y, while if a^-i = — 1, 4rfe +1 is within | of 
|, so that < 2rfc + i < ^ + |, and | — ^ < r^+i < i. In the first case 
we will make use of the fact that as y decreases by a factor of |, S decreases 
by a factor of In the second case we make use of the factor hij(z) with 
|| £ — A || < 4, where || • || denotes the distance from the nearest integer. We 
note that in this second case, hij(z) is necessarily h\ t \{z) since the other filters 
are at such a z, and we were assuming that we were in the support of 4>i(x). 
Thus, in the second case, we obtain a factor of \p{z)\ with \\z — ^|| < j, and 
we reset 8 to jg. Since x e A n U B n U C„ satisfies 4j- < \x\ < 8" +1 , it will 
take between n and n + 1 applications of this process of moving from y to | 
to move into the unit interval. From this, if aj = 1 for all j, we will obtain 
n factors of the form p{z) for ||z — A|| < j^. Thus, by applying the Mean 
Value Theorem to condition (5) in the definition of the function p, we see that 
we will have n factors |/ii,j(fr)| each of which is less than (jj2") r+2 Each time 
a 



= +1, wc will decrease the value of 5 by a factor of | instead of obtaining 
a new small factor. If aj = +1 for all j, we will eventually arrive at ^1^(2;) 
for |^ — i| < ■^ r . Again, by applying the Mean Value Theorem to condition 
(5) in the definition of p, we see that this gives us a factor of |^jj(^-)| which 
is less than (gt) r+2 - Values of x G A n U B n U C n with a mixture of aj = +1 
and a j ■ — — 1 will be somewhere between these two extremes, but all will have 
one or more factors of p(x) among the |/iij(|r)l with x so close to or j that 
the product of these factors is less than (j^r) r+2 . Since the other factors are 
bounded by 1, we thus get that each term in 4*i(x) is bounded in absolute value 
by 7f 3n (^r)'' +2 - Thus, for x E A n U B n U C n , we have \4>i(x)\ < ^f^^ ■ 

Finally, for such x, \x r+1 <fii(x)\ < ^^jf^+q'^ = 3 " +1 ^'' +1 — > as n — > 00. □ 
We are now ready to prove our main theorem: 

Theorem 4. For any fixed integer r > 0, there exists a real-valued normalized 
tight frame wavelet i\) such that tp € C r , ip € C°° and such that the multiplicity 
function associated to ip is that of the Journe wavelet. 

Proof. We build tf> using the technique described in Theorem |3 with low-pass 
filter functions hi j defined in equations (|13[) - l|16|l . and with high-pass filters 
given by 



9i{x) 



92{x) 



e 2 "p(i + |) xe[-f,f) 
otherwise 



~e 2mx p{x) x<E[=±,\) 
otherwise 



It is easily checked that these generalized high-pass filter functions satisfy 
the filter equations (|llfl and l|12l) . In fact, as their form suggests, they were built 



1G 



from the generalized low-pass filters using the classical high-pass filter obtained 
from p via equation . The frame wavelet ip is defined as in Theorem |3J by 

ip(x) = (ffi(§)0i(§) + 32(f )02(§ j where 0i and 02 are the functions built 
from the {hi j} in Proposition^ The wavelet ip is real- valued since the fact that 
both the low and high pass filters are symmetric forces ip to be symmetric as 
well. We have that ip is C°° since Lemma[3]established that 0i and 02 are C°°, 
and gi and g 2 are as well. Further, we have that |a:| r+1 '0(a;) — > as x — > 00, 
since 0i and 02 were shown to have this property in Proposition Eland g\ and 
g 2 are bounded by y/2. Thus the wavelet "0 itself satisfies tp € C r . 

It remains to show that ip determines the same GMRA as 0i and 02, so 
that this wavelet also has the multiplicity function of Journe. It was shown 
in [H] that ip <E Vi, where {Vj} is the GMRA determined by 0i and 02 as in 
Proposition^] This says that ip is obtained from the GMRA in the sense defined 
by Zalik 2^j. We will now show that in this particular case, ip is associated with 
the GRMA {Vj} in the sense that Vj is the closed linear span of {ipi,k}i<j- (For 
more on this distinction, see e.g. JOj-) Since we already know (by ip £ Vi) that 
Vj contains this closed linear span, it will suffice to show that Vq is contained 
in the closure of the span of {ipi,k}i<o- 

We will do this by defining a bounded surjective linear map T : L 2 {S\) © 
L 2 {S 2 ) i-> V , that takes a tight frame for L 2 (Si) © L 2 {S 2 ) to {ipi,k}i<a- We 
define T(/i, f 2 ) — /i0i + f 2 (p2, where fi £ L 2 (Si) is extended to a periodic func- 
tion. T maps onto Vq since by Lemma ^ the support of (pi is the periodization 
of Si. To see that T is bounded, we compute the norm 



\\n{h})\\ 2 = / hWMx) + h( X )<h{x) 



dx 



fi{x)(pi(x) + f 2 (x)<p 2 (x) 
E [ \h(x)(pi(x + j) + f 2 (x)(p 2 (x + j)\ dx 

/ i/i(s)i a E \m* + j)f + E \m x + 3) 



dx, 



where the second equality follows from the disjointness of the supports of 0i 

^2 

and 02- Let Per0j(x) = Sjez 4>i{ x + j) ■ To show that T is bounded, it will 

suffice to show that Per0^ is bounded on T for i = 1,2. This follows from 
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Lemma n by writing, for x € [— h, i] 



Per^z) = ^ M x + ^2 a j 8j ) 

„=0{ aj }e{-i,i}" 3=0 



Z-^i g(n— l)r 
n— 

where the last step follows from Proposition 

Now to see that T takes a tight frame for L 2 (Si) © L 2 {S 2 ) to {^;,fc}/<o, we 
adopt some notation from 0. Define the operator Sg '■ L 2 (T) — > L 2 {S\)®L 2 {S2) 

by 

[S G (f)](x) = { 9l (x)f(2x),g 2 (x)f(2x)} 
and 5 H : L 2 (S*i) © L 2 (S* 2 ) -» L 2 (Si) © L 2 (5 2 ) by 

[5 ff ({/i,/ 2 })](x) = {ftiA(s)/i(^)+'»2 > i(a:)/3(2a:),^(a;)/i(2a;)}. 

It is shown in (Lemma 3.3) that {S%{S G (e 27rikx ))} keZtn=0A , 2 ,.. forms a nor- 
malized tight frame for L 2 (Si) © L 2 {S 2 ). We complete the proof by showing 
that the operator T takes this frame to {ipi,k}i<o- Note first that for n = 0, wc 
have 

T{S G (e 2mkx )) = 9l (x)e^ kx Mx) + g 2 (x)e^ kx Mx) 
= e 4nkx $(2x) 
= i>-l,k 

Thus the result will follow from the fact that T intertwines the operators Sh 
on L 2 (Si) © L 2 (S 2 ) and inverse dilation S^ 1 on V : 



[T(Sh({A, h}))](x) = h 1A (x)h(2x) +Ji 2A (x)f 2 (2x)Mx) + hi fl {x)fx{1x)ihi^ 

= ,h(2x) (h lt i(x)Mx) + hiA x )M%)) + h{2x)h2,i{x)4> 1 {x) 
= h{2x)V2$i{2x) + f 2 (2x)V2fo{2x) 



S- 1 (h4>i + h4> 2 

5-\T{f))]{x). 



(x) 



□ 



Remark 3. While the dilates and translates of the wavelet ip defined in Theorem 
^] do form a normalized tight frame for all of L 2 (R), the translates of ip do not 
even form a frame for their linear span. We can see this using the result in 



once again, since Per -0(2;) = ^ 



is not bounded away from on 



ip{x+j) 

its support. In particular, Per ip is continuous, yet nonzero on ±(f, f), and 
on±(f,f). 



18 



References 



[1] A. Aldroubi, C. Cabrclli, and U. Molter, Wavelets on irregular grids with 
arbitrary dilation matrices and frame atoms for L 2 (M. d ), Appl. Comput. 
Harmon. Anal. 17, 119-140 (2004). 

[2] P. Auscher, Solution of two problems on wavelets, J. Geom. Anal. 5, 181— 
236 (1995). 

[3] L. W. Baggett, Redundancy in the frequency domain, Harmonic Analysis 
and Applications (C. Heil, ed.), Birkhauser, Boston, to appear. 

[4] L. W. Baggett, J. E. Courter, and K. D. Merrill, The construction of 
wavelets from generalized conjugate mirror filters in L 2 (W l ), Appl. Comput. 
Harmon. Anal. 13, 201-223 (2002). 

[5] L. W. Baggett, P. E. T. Jorgensen, K. D. Merrill, and J. A. Packer, Con- 
struction of Parseval wavelets from redundant filter systems, J. Math. Phys. 
46 (2005), no. 8, 083502, 28 pp. 

[6] L. W. Baggett, H. A. Medina, and K. D. Merrill, Generalized multi- 
resolution analyses and a construction procedure for all wavelet sets in 
R™, J. Fourier Anal. Appl. 5, 563-573 (1999). 

[7] L. W. Baggett and K. D. Merrill, Abstract harmonic analysis and wavelets 
in R™, in The Functional and Harmonic Analysis of Wavelets and Frames 
(San Antonio, TX) (L. Baggett and D. Larson, eds.), Contemp. Math. 247, 
Amer. Math. Soc, Providence, RI, 1999, pp. 17-27. 

[8] J. Benedetto and M. Leon, The construction of multiple dyadic minimally 
supported frequency wavelets on R d , in The Functional and Harmonic 
Analysis of Wavelets and Frames (San Antonio, TX) (L. Baggett and D. 
Larson, eds.), Contemp. Math. 247, Amer. Math. Soc, Providence, RI, 
1999, pp. 43-74. 

[9] J. J. Benedetto and S. Li, The theory of multiresolution analysis frames 
and applications to filter banks, Appl. Comput. Harmon. Anal. 5, 389-427 
(1998). 

[10] M. Bownik, Riesz wavelets and generalized multiresolution analyses, Appl. 
Comput. Harmon. Anal. 14, 181-194 (2003). 

[11] O. Bratteli and P. Jorgensen, Wavelets Through a Looking Glass: The 
World of the Spectrum, Birkauser, Boston-Basel-Bcrlin, 2002. 

[12] A. Cohen, Wavelets and Multiscale Signal Processing, translated by R. 
Ryan, Chapman and Hall, U.K., 1995. 



19 



[13] J. Courter, Construction of dilation d wavelets, in The Functional and 
Harmonic Analysis of Wavelets and Frames (San Antonio, TX) (L. Baggett 
and D. Larson, eds.), Contemp. Math. 247, Amer. Math. Soc., Providence, 
RI, 1999, pp.183-205. 

[14] J. Cuntz, Simple C*-algebras generated by isometries, Comm. Math. Phys. 
57, 173-185 (1977). 

[15] X. Dai, D. R. Larson, and D. Speegle, Wavelet Sets in W 1 , J. Fourier Anal. 
Appl. 3, 451-456 (1997). 

[16] I. Daubechies, Ten Lectures on Wavelets, American Mathematical Society, 
Providence, RI, 1992. 

[17] I. Daubechies, Orthonormal bases of compactly supported wavelets, Comm. 
Pure Appl. Math. 41, 909-996 (1988). 

[18] A. Haar, Zur theorie der orthogonalen funktionen systems, Math. Ann. 69, 
331-371 (1910). 

[19] D. Han, D. R. Larson, M. Papadakis, and Th. Stavropoulos, Multircsolution 
analyses of abstract Hilbert spaces and wandering subspaces, in The Func- 
tional and Harmonic Analysis of Wavelets and Frames (San Antonio, TX) 
(L. Baggett and D. Larson, eds.), Contemp. Math. 247, Amer. Math. Soc, 
Providence, RI, 1999, pp. 259-284. 

[20] E. Hernandez and G. Weiss, A First Course on Wavelets, CRC Press, Boca 
Raton-New York, 1996. 

[21] W. M. Lawton, Tight frames of compactly supported affine wavelets, J. 
Math. Phys. 31, 1898-1901 (1990). 

[22] S. Mallat, Multiresolution approximations and wavelet orthonormal bases 
of L 2 (R), Trans. Amer. Math. Soc. 315, 69-87 (1989). 

[23] Y. Meyer, Wavelets and Operators, Cambridge Studies in Advanced Math- 
ematics v. 37, Cambridge University Press, Cambridge, England, 1992. 

[24] R. Strichartz, How to make wavelets, Amer. Math. Monthly 100, 539-556 
(1993). 

[25] E. Weber, The action of translations on wavelet subspaces, Ph.D. Thesis, 
University of Colorado (1999). 

[26] R. A. Zalik, Riesz bases and multiresolution analyses, Appl. Comput. Har- 
mon. Anal. 7, 315-331 (1999). 



20 



